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ABSTRACT 

We estimate the size and distribution of the parent populations for the 6 largest (at 
least 20 stars in the Solar neighborhood) chemical groups identified in the Chemical 
Tagging experiment by Mitschang et al. 2014. Stars in the abundance groups tend to 
lie near a boundary in angular momentum versus eccentricity space where the proba¬ 
bility is highest for a star to be found in the Solar neighborhood and where orbits have 
apocenter approximately equal to the Sun’s galactocentric radius. Assuming that the 
parent populations are uniformly distributed at all azimuthal angles in the Galaxy, 
we estimate that the parent populations of these abundance groups contain at least 
200,000 members. The spread in angular momentum of the groups implies that the 
assumption of a uniform azimuthal distribution only fails for the two youngest groups 
and only for the highest angular momentum stars in them. The parent populations of 
three thin disk groups have narrow angular momentum distributions, but tails in the 
eccentricity and angular momentum distributions suggest that only a small fraction 
of stars have migrated and increased in eccentricity. In contrast, the parent popula¬ 
tions of the thick disk groups exhibit both wide angular momentum and eccentricity 
distributions implying that both heating and radial migration has taken place. 


1 INTRODUCTION 


|Freeman fe Bland- Hawthorn (2002) proposed that stars with 
similar abundance measurements could represent a particu¬ 
lar star formation or enrichment event, discrete in space and 
time, as might be expected from the homogeneity of nearby 
open clusters and moving groups ( De Silva et al.|2006 


2007 


Quillen 2002). These stars would subsequently disperse 


the Galaxy, retaining their initial chemical patterns (Bland- 
|Hawthorn et al.||2010 ). A search for stars that have a sim¬ 
ilar abundance pattern as the Sun would allow us to learn 
about the birth place of the Sun (|Portegies Zwart 2009) |Liu 
et al.|2015 ). Because stars in the solar neighborhood span a 


wide distribution in stellar ages, metallicities and inferred 
birth Galactocentric radii, it is difficult to pin down the 
role of specific mechanisms for stellar migration and heat¬ 
ing (increase in radial and vertical epicyclic motions) (e.g., 


see|Quillen et al.||2009 Schonrich &; Binney||2009| |Fr> 


&; Bland-Hawthorn|2002[ |Kruijssen et al.|2011| Haywood et 


al.||2013| |Minchev et al.|2014| Lehnert et al.|2014| ). A study 

of homogeneous groups of stars should give complementary 
constraints on migration and heating processes, compared 


to those arising from studies of heterogeneous distributions 
(such as a magnitude limited sample of stars in the Solar 
neighborhood). 

We focus on the a high-resolution spectroscopic study 
of 714 F and G dwarf and subgiant stars in the Solar neigh¬ 
borhood studied by Bensby et al. (2014). The blind chemi¬ 


cal tagging experiment by Mitschang et al. (2014) used this 


sample to identify groupings of nearby disc field stars that 
share metal abundance measurements. The field stars they 
identified as having similar abundances are not clustered 
in space, nor do they share similar space motions. Using 
isochrone sets, Mitschang et al. (2014) estimated the ages of 


each of these chemical groupings. These groups represent a 
first attempt to identify groups of stars from single discrete 
birth events. 

We ask here: what is the number and distribution in the 
Galaxy of a parent stellar population of one of these abun¬ 
dance groups? We necessarily focus on only the 6 largest 


groups identified by Mitschang et al. (2014) each of which 


contains more than 20 stars. We begin by assuming that 
the parent population for each group is currently evenly 
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2 Quillen et al. 


distributed (azimuthally) in the Galaxy and at the cur¬ 
rent time only a fraction of the stars in the parent popu¬ 
lation are present in the Solar neighborhood. This assump¬ 
tion was adopted for the toroid models by Bland-Hawthorn 
et al. ( 2010| ) (for an illustration see their Figure 4). This as¬ 
sumption neglects how a cluster dissolves and is dispersed in 
the Galaxy (see discussions by Bland-Hawthorn et al.| 2010 


Portegies Zwart|2009 ). We will discuss how our assumption 
of axi-symmetry for the parent distribution could have im¬ 
pacted our inferred parent population distributions. We also 
neglect our location in the Galaxy with respect to spiral and 
bar perturbations, (e.g., Quillen et al. 2011 Quillen 2014 


Minchev et al. 2013, 2014). We assume there is no correla¬ 


tion between vertical oscillation amplitude and eccentricity 
in the parent population and that the epicyclic angle distri¬ 
bution is relaxed (see Minchev et al. 2009] for an illustration 
of what can be seen when this is not true). 

We first consider the fraction of time that a star with 
a given eccentricity and angular momentum might be seen 
in the Solar neighborhood [j From a distribution of orbits 
with a given eccentricity and angular momentum we esti¬ 
mate the probability that a star is seen in a Solar neigh¬ 
borhood volume with boundary dependent on the distances 


of stars in the Bensby et al. (2014) sample. For each star 


in one of Mitschang et al.’s chemical groupings the inverse 
of this probability lets us estimate the number of stars at 
similar eccentricities and angular momentum in the parent 
population. 


Bensby et al. (2014) selected stars for spectroscopic 


study with a range of properties and necessarily did not 
observe every F and G star in the Solar neighborhood. We 
compare the Bensby et al. (2014) sample with the Geneva- 


Copenhagen survey of F and G stars in the Solar neighbor¬ 


hood (Nordstrom et al. 2004 Holmberg et al.|[2009 ) to es¬ 
timate a selection bias as a function of angular momentum. 
The GCS is a magnitude limited, kinematically unbiased 
sample of almost 17000 nearby F and G stars. [Bensby et ah] 
(2014) warn that their sample is a compilation of a number 
of different observing programs and so they give no selection 
description for the entire sample. We should be careful inter¬ 
preting inferred parent distributions, keeping in mind that 
there might be additional biases arising from the selection 
of this sample. 

Using these two corrections, the first based on probabil¬ 
ity for such an orbit to be seen in the Solar neighborhood, 
the second based on selection bias, we derive estimates for 
the source or parent populations of the 6 abundance groups 
identified by Mitschang et al. (2014). A discussion follows 


on the nature of the parent populations and on how our un¬ 
derlying assumptions have impacted our estimate of their 
number and distributions. 


1 We only consider the z component of angular momentum that 
is dominated by rotation in the Galaxy. In the Solar neighborhood 
the angular momentum of a star L « Rq(V + Vlsr,q) with Rq 
the galactocentric radius of the Sun, Vlsr the rotation velocity 
of the local standard of rest and V the tangential component of 
the star’s velocity vector. 


2 PROPERTIES OF THE SIX ABUNDANCE 
GROUPS 

Properties of the 6 largest abundance groups found by 
Mitschang et al. (2014) are listed in Table [l] The Hippar- 


cos catalog numbers of the stars in each group are listed in 
Table 2 by Mitschang et al. (2014). Our Table [l] lists the 


group identification number (from their Table 2), group age 


(that derived by Mitschan g et al.| 2014 using Yonsei-Yale 
isochrone sets and with error estimate described in their 
section 4.2) the mean [Fe/H] of the group (in Solar units 
and using abundances listed in Table C2 by |Bensby et al.| 


2014) and the group mean [a/Fe] in Solar units. The mean 


abundance values and standard deviation for each group are 
computed from the values for each star in the group. For 
each star [a/Fe] is calculated by averaging the abundances 


for a elements Ti, Mg, Si, and Ca (as done by Mitschang 
et al. Table a also lists the mean (L) and standard 


deviation, ol, (in km s -1 kpc) of the angular momentum 
distributions of each group. These are the angular momen¬ 
tum values by B ensby et al.| (2014) who computed space 
motions for all the stars in their sample (see their section 

3)Q 

Studies of abundance populations based on high resolu¬ 
tion spectroscopy find a bi-modality in the abundance distri¬ 
bution (e.g.,|Navarro et al.|2011 


Fuhrmann 2011 [ |Haywood 
|et al.|2013||Adibekyan et al.|2013|| Anders et al.|2014| and ref- 

erences therein) with a dividing line between thin and thick 
disks populations near [a/Fe]~ 0.12 (e.g., see Figure 12 by 
Reddy et al.|2006 and section 5.1 by |Mitschang et al.|2014 ). 


For each abundance group, mean values of [Fe/H] are plotted 
against the mean values of [a/Fe] in Figure [l] with the other 
stars in the Bensby et al. ( 2014| sample. The three youngest 
groups have abundances consistent with a thin disk popu¬ 
lation, whereas the older three have abundances consistent 
with a thick disk population. 


690 out of 714 stars in the Bensby et al. (2014) sample 


are also present in the Geneva-Copenhagen Survey of F and 
G stars in the solar neighborhood (GCS) (Nordstrom et al 


2004 Holmberg et al. 2009). The distribution of distances 


from the Ben sby et al.|( 2014) sample is compared to that 
of the GCS in Figure [2 illustrating that the stars in the 


Bensby et al. (2014) sample are predominantly nearer than 


100 pc. Here distances are based on parallaxes from the new 
reduction of the Hipparcos data by van Leeuwen (2007). Of 


the 163 stars in the 6 abundance groups, we find that only 
9 of the stars are further than 100 pc from the Sun. Thus 
this sample of stars is confined to a small spherical volume, 
centered on the Sun, with an approximate radius of 100 pc. 

From the angular momentum, L, and eccentricity, e, 
values listed by Bensb y et al.| ( |2014 ) we constructed his¬ 
tograms for each group, and these are shown in Figure[3] The 
mean angular momentum (also listed in Table |T]) for each 
group decreases with increasing age, suggesting that the old¬ 
est groups arise from the inner galaxy and the youngest 
groups are located near the Sun’s galactocentric radius. 


While Bensby et al. (2014) did not list errors for ec¬ 


centricity e or angular momentum L for each star, we can 


2 The adopted local standar d of rest ((7 q, Vq, VFq) = (11-10, 
12.24, 7.25) km s -1 is that by Schonrich et al. (2010). 
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Parent populations of abundance groups 3 



Figure 1. Mean abundances (in solar units) for each of the six abundance groups (found by |Mitschang et al.~1 |2014| are shown as 
large red dots. Error bars represent standard deviations of the abundance values of stars in the group. Black dots show abundances of 
individual stars from the Bensby et al.| (2014| sample. The dashed lines show solar values. The three youngest groups are typical of the 
thin disk population, whereas the three oldest are typical of the thick disk population. 
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Figure 2. Distance distribution for thejBensby et al. (2014|) sample compared to that of the Geneva-Copenhagen Survey ([Nordstrom et| 
|al.|2004||Holmberg et ah| 2009). Most stars from the [Bensby e t al. (2014J sample are within 100 pc of the Sun. The distributions have 
been normalized so that they integrate to 1. 


assume that the space velocity components U, V, W have er¬ 
rors the same size as those of the GCS wh ich are estimated 
to be Av ~ 1.5 km s -1 (see section 4.7 by Nordstrom et al. 


2004). This corresponds to an approximate error of AL ~ 13 
km s _1 kpc in angular momentum. We estimate the size of 
an error in eccentricity with Ae ~ A v/Vlsr,q ~ 0.01. The 
errors could also have systematic trends in them (as a func¬ 
tion of other parameters such as position on the sky) and 


due to uncertainty in the Solar motion or the rotation curve 
used to calculate the eccentricity^] 


3 At low eccentricity the difference E — E(L) ~ 0.5 ft 2 r 2 e 2 
where r g is the guiding radius, e is the orbital eccentricity, E 
is the orbital energy per unit mass, E(L) is the energy (per unit 
mass) of a circular orbit with angular momentum L, and ft is the 
epicyclic frequency. For a power law rotation curve v c (r ) oc r~ a 
the epicyclic freqency ft = -\/2(l — a)Q where the angular rota- 
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4 Quillen et al. 


Table 1 . Properties of the Abundance Groups 


GID 

N 

Age 

(Gyr) 

[Fe/H] 

[a/Fe] 

(L) 

(km s 

°L 

1 kpc) 

<e> 

5 

21 

4.0 ±0.6 

0.23 ± 0.04 

0.01 ±0.02 

1715 

214 

0.14 

1 

42 

4.8 ±0.4 

0.07 ±0.07 

0.01 ±0.01 

1748 

262 

0.17 

3 

25 

7.1 ±0.4 

-0.07 ±0.06 

0.03 ± 0.02 

1522 

247 

0.28 

4 

24 

10.1 ± 1.4 

-0.43 ± 0.05 

0.18 ±0.02 

1450 

271 

0.32 

2 

30 

10.2 ±0.8 

-0.30 ±0.05 

0.22 ±0.02 

1368 

373 

0.34 

6 

21 

12.1 ± 1.1 

-0.64 ±0.07 

0.26 ± 0.02 

1297 

510 

0.39 


GIC is the group number given by Mitschang et al.| ( |2014| >. N is the number 
of stars in the abundance group and the estimated age by [Mitschang ^et al.| 
§014^ is given in Gyr. Mean [Fe/H] and [a/Fe] values in solar units for the 
stars in the group were computed using abundance values listed in Table C.3 
by |Bensby et ah] ( |2014| >, finding them using Hipparcos catalog numbers for 
the stars listed in Table 2 by [Mitschang "et al.| ( |2014| for each group. Errors 
in the abundances are the standard deviations of the abundance values from 
each star. We list the mean (L) and standard deviation, a l, (in km s —1 kpc) 
of the angular momentum distributions and the mean eccentricity (e) for 
each group (distributions are shown in Figure |3|. These are computed from 
the eccentricity and angular momentum values computed by |Bensby et al.| 
§014). 


3 PROBABILITY OF DETECTING A STAR IN 
THE SOLAR NEIGHBORHOOD AS A 
FUNCTION OF ORBITAL ECCENTRICITY 
AND ANGULAR MOMENTUM 


In this section we describe how to estimate the probability 
that a star with a given angular momentum and eccentricity 
is found in a solar neighborhood sample if the parent pop¬ 
ulation is randomly distributed in azimuthal and epicyclic 
angles. We use angular momentum and eccentricity to de¬ 
scribe each orbit. We use eccentricity e instead of energy 
as it is unit less, describes the extent of radial excursion in 
the orbit and so gives an intuitive description for the orbit 
shape, and it does not depend on a potential energy offset. 

To be consistent with the angular momentum and ec¬ 


centricities computed by Bensby et al. (2014) and Mitschang 


et al. (2014) we use the same Galactic potential model as 


they did to compute our probabilities. Using the gravita¬ 
tional potential for the Galactic model by Allen & Santillan 
(199lf we integrate planar orbits with different initial an¬ 


gular momentum and different initial radii. For each orbit we 
record the eccentricity defined as e = (R a — R p )/(R a ± R P ) 
(following |Bedin et al.| 2006) where R a ,R P are radii of galac¬ 
tic apocenter and pericenter, respectively. 

For each e, L, we computed a few thousand positions 
in a full orbit (using a finite size timestep to compute a 
full orbital period). We then randomly chose a few thou¬ 
sand azimuthal angles (corresponding to randomly chosen 
initial orientations) giving a total of approximately 10 mil¬ 
lion points in the galaxy plane to compute each probability. 
At each timestep and for each angle we computed the po¬ 


tion rate Q = v c /r. Uncertainty in the slope of the rotation curve 
affects the estimate for E(L) and epicyclic frequency, r, and so 
the computed values for the eccentricity. 

4 This model assumes a Galactocentric distance for the Sun and 
rotation velocity of a circular orbit at that radius of Rq =8.5 kpc 
and Vlsr,q = 220 km s -1 . With these values the angular mo¬ 
mentum of the local standard of rest is Llsr = 1870 km s -1 kpc. 


sition of the star and the fraction that fell within a solar 
neighborhood area, within 100 pc of the Sun, gave the prob¬ 
ability. For randomly distributed initial azimuthal angle and 
initial position in the orbit, and using a single orbital period, 
we measure the fraction of stars in an orbit, as a function of 
angular momentum and eccentricity, that are located within 
100 pc of the Sun. In other words, we assume there is a dis¬ 
tribution of orbits with this angular momentum and eccen¬ 
tricity that is randomly distributed in azimuthal angle, and 
using this distribution, we compute the probability, p 0 (e, L), 
that a star would be observed in the Solar neighborhood at 
any particular time. For a range of orbital eccentricities, e, 
and angular momenta, L, we compute p 0 (e, L) and display 
it in Figure [3] The color bar shows the log 10 of the prob¬ 
ability. The black dots show the orbits that we integrated 
and that were used to make the color contours. The angular 
momentum is in units of km s -1 kpc. Wiggles in Figure [3] 
are artifacts due to the sampling of the orbits integrated. 

At low eccentricity and angular momentum L above or 
below that of the local standard of rest, the probability p 0 , 
for the orbit is zero as the orbit never crosses the Sun’s 
galactocentric radius. A star in such an orbit is never near 
the Sun. The white region on the lower left and upper left 
in Figure is this forbidden region. Large eccentricity or¬ 
bits that do cross the Sun’s galactocentric radius (on the 
right in Figure [3]) are less probable than lower eccentricity 
ones as stars spend much of the time at larger or smaller 
galactocentric radius than that of the Sun. For a given an¬ 
gular momentum, the probability is highest at an eccentric¬ 
ity that just barely allows the orbit to cross into the Solar 
neighborhood. We attribute the increase in probability near 
the forbidden region boundary to the large fraction of the 
orbital period spent near a particular radius when at apoc¬ 
enter or pericenter. This effect has previously been described 
as a bias due to crossing times in the Solar neighborhood 
). The effect is illustrated in Figure [3] 
showing epicyclic oscillations for three different groups of 
orbits, one with apocenter near the Sun’s galactocentric ra¬ 
dius that is likely to be seen in the Solar neighborhood, high 


(Mayor et al. 1977 
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Figure 3. Eccentricity and angular momentum distributions of the 6 groups. The left panels show the eccentricity distributions and the 
right panels show the angular momentum distributions. Each row shows a different group, with group labelled by age. Shown in blue 
with axis on the right hand side of the left panels are the observed eccentricity distribution of each group, plotting numbers of stars 
in bins of width 0.05 in eccentricity. Shown in green, with axis on the right hand size of the right panels are the angular momentum 
distributions of each group, showing the number of stars in bins of width 100 km s —1 kpc. The pink shaded regions, corresponding to 
axes on the left side of each panel, show the estimated numbers stars in the parent population distributions (corrected for selection and 
assumed to be evenly distributed in azimuthal angle, see Section 4) with bin widths of 0.01 in eccentricity or 20 km s —1 kpc in angular 
momentum. 


eccentricity orbits that have a lower probability and orbits 
within the forbidden region that cannot be found in Solar 
neighborhood. 

We can account for this probability increase near apoc- 
enter using an epicyclic approximation for radial orbital 
variations. For low eccentricity stars the radius r{t) ~ 
r g (l + ecos(ft£ + 0o)) where r g is the guiding radius, e is the 
eccentricity, k is the epyclic frequency, (f>o an initial phase 
and the apocentre radius R a — r g ( 1+e). Near apocenter and 
using a small angle approximation, R a — r(t) oc (t — t apo ) 2 
where t apo is a time when the orbit is at apocenter. This 


gives a dependence of the fraction of the orbital period, /, 
spent within a narrow annulus of width dr from apocenter, 
/ oc y/d/r. In contrast when the orbit is near the guiding ra¬ 
dius and using a small angle approximation, r — r g oc {t — t g ) 
(with t g a time the orbit crosses the guiding radius) giving a 
dependence of the fraction of the orbital period spent within 
dr of r g to be / oc dr. For a small range of radius dr , the 
fraction of the orbital period spent near apocenter is larger 
than that spent near the middle of the orbit at the guiding 
radius. The trend is still present at moderate eccentricity 
where the epicyclic approximation is less accurate. 
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Figure 3. continued 


3.1 Probabilities of individual abundance group 
stars 

We now consider the probabilities that stars in the abun¬ 
dance groups are seen in the Solar neighborhood. In Figure 
[6] we show the eccentricity and angular momentum of stars 
in these abundance groups on top of the probability, p 0 , that 
such a star is found in the Solar neighborhood. Each panel 
shows a different abundance group and the groups are la¬ 
belled by their ages. The probability is displayed as in Figure 
[4] Figure [6] illustrates that stars tend to be found near the 
forbidden region in e, L space, as expected from the location 
of the high values in probability distribution p G (e, L). 

For each star with eccentricity e* and angular momen¬ 
tum Li we can use the probability p G (ei, Li) to estimate the 
size of the parent grouping. The parent population has at 
least 


*i = £ 

i 


1 

Po(ei, Li) 


( 1 ) 


stars in it. If we underestimate the probability p 0 (ei,Li) 
then we will overestimate the number of stars in the parent 
population. To ensure that observational errors in a and 
Li for individual stars do not give spurious high numbers 
near the forbidden region, we take p Q to be the maximum 
value within e* d= Ae and Li ± A L with Ae = 0.01 and 
AL = 13 km s -1 kpc, the size of the errors estimated for 
these quantities (see end of section 2). For each abundance 
group, we have summed the inverse of the probabilities and 
list the total number of estimated parent stars in Table [2] 
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Figure 4. As a function of orbital eccentricity and angular momentum we show log 10 p o (e,L) of the probability that a distribution 
of stars in such an orbit (but with randomly chosen angles) would be observed in the Solar neighborhood (100 pc from the Sun). The 
probability is computed using planar orbits integrated in the Galactic potential model by [Allen & Santillan (1991). Each orbit integrated 
is shown as a black dot. The angular momentum is in units of km s —1 kpc. 
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High eccentricity orbits crossing the Sun's galactocentric radius 



Figure 5. Orbits with apocenter near Sun’s galactocentric radius (as shown on the top) are more likely to be seen in a solar neighbor¬ 
hood sample, than a high eccentricity orbit (as shown in the middle). The bottom illustrates orbits that never are found in the Solar 
neighborhood, corresponding to a region in e, L space that we call the forbidden region. 
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Figure 6. Eccentricity and angular momentum of stars in an abundance group plotted as points on top of the probability that such an 
orbit is detected in the Solar neighborhood. Each panel shows a different abundance group and the groups are labelled by their ages. 
The color bar shows log 10 p o of the probability computed for each orbit. The stars in each abundance group often he in regions of high 
probability, near the boundary of the forbidden region. 


4 PARENT POPULATION DISTRIBUTIONS 

Using the probabilities estimated for each star in an abun¬ 
dance grouping we now estimate the number density of stars 
in the grouping as a function of eccentricity and angular mo¬ 
mentum e, L. Each star with e*, Li contributes a total parent 
population of p 0 (ei, Li) -1 at e^L*. We smooth this distri¬ 


bution to estimate the number density of stars in the parent 
population as a function of e, L. For each cluster, the re¬ 
sulting distributions in e and L are shown in Figure [7] The 
color bars show the number of stars per eccentricity and an¬ 
gular momentum bin with bin size de = 0.01 and dL = 20 
km s _1 kpc. The distributions have been smoothed by 4 or 
5 eccentricity bins and 2 or 3 angular momentum bins with 
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the tighter distributions (for the younger groups) smoothed 
by fewer bin widths. 

Bensby et al. ( 2014| selected stars for spectroscopic 


study so that the sample contained extremes of both thin 
and thick disk. Therefore, many thin disk stars were neces¬ 


sarily neglected from the Bensby et al. (2014) sample. As 


a result, the sample contains a bias against high angular 
momentum stars. In Figure |8| we compare the angular mo- 


1L 

;|Be 


mentum distribution o f the|Bensby et al. (2014) sample to 
that of the GCS stars (Hol mberg et al.|2009 ), but restricted 
to stars within 80 pc. In this figure the blue histogram shows 


the Bensby et al. (2014) sample and the red histogram the 


GCS stars with overlaps displayed as purple. 

We constructed a selection function, /(L), choosing two 
tanh functions as they go smoothly between one constant 
to another constant value and this allows us to model the 


the two humps in the Bensby et al. (2014) sample evident 


in the angular momentum distribution shown in Figure [5] 
The function we chose is described with a few parameters, 
is smooth, is never extremely small and does not cross zero 
(this is important as we need to divide by it). The black dots 
in Figure[8]show the GCS histogram multiplied by /(L) with 

- 7 — ) - tanh ( - 7 — ) ( 2 ) 

siLlsr J 2 \S2-LlSR J 


f(L) = a 0 - y tanh 


and coefficients do = 0.53, a\ — 0.65, — 0.30, si = 

0.18, S 2 = 0.05 and L\ — 1148, L 2 = 1445 km s -1 kpc and 
with Lslr = 1870 km s -1 kpc. The function itself is plotted 
in Figure[9] We did not automatically fit the coefficients, but 
did adjust the coefficients so that the two histograms lay on 
top of one another. Because we divide by this function, if 
it is an underestimate for the selection of the Bensby et al. 
sample, then we will overestimate the size of parent popula¬ 
tions. 

The function f(L) is an estimate for the fraction of 
stars selected by |Bensby et al. (2014) compared to that in 
the Geneva Copenhagen Survey. At L near that of the lo¬ 
cal standard of rest, f(L) ~ 0.05 implying that for every 


high L star in the Bensby et al. (2014) sample there are 20 
stars with similar angular momentum in the GCS. At high 


angular momentum, the Bensby et al. (2014) sample also 


contains a higher proportion of high eccentricity stars than 
the GCS. To ensure that we don’t overestimate the num¬ 
ber of high angular momentum high eccentricity stars the 
parent populations we cut the selection function with 

/(L,e) = 1 for L > 1600 km s -1 kpc and e > 0.15 


= /(L) otherwise 


(3) 


We correct the probability for each star in each abun¬ 
dance grouping with this selection function giving a total 
number for the parent population 


i 


_1_ 

p 0 (e-i : Li)f{Li , ei) 


(4) 


again taking the minimum value for p Q within e* d= Ae and 
Li =b AL. The distribution of the parent populations for 
each abundance grouping, also taking into account the se¬ 
lection function, are shown in Figure [To] and the numbers 
N 2 listed in Table [2] As expected the total number of stars 
estimated for each group is larger than that estimated pre¬ 
viously without using the selection function. We integrate 
the parent distributions plotted in Figure in eccentric¬ 


ity to estimate angular momentum distributions. Likewise 
integrating in angular momentum we can estimate eccen¬ 
tricity distributions. The estimated parent eccentricity dis¬ 
tributions and the parent angular momentum distributions 
are shown as pink solid regions in Figure [3] where numbers 
of stars in eccentricity bins of size 0.01 are plotted in the 
left panels and numbers of stars in angular momentum bins 
of size 20 km s -1 kpc are plotted in the right panels. 


4.1 Discussion on azimuthal structures and phase 
wrapping 

To estimate the probability p a we have assumed that the 
parent population is evenly distributed azimuthally in the 
galaxy. However, an originally cold disrupted cluster may 
not have time to become evenly distributed in azimuthal an¬ 
gle (for illustration see Figure 2 by Po rtegies Zwart| |2009). 
This would lead to a bias - a survey of the Solar neighbor¬ 
hood would not see every group that is present at the Sun’s 
galactocentric radius. Furthermore, the parent population 
size of a group detected in the Solar neighborhood would be 
overestimated by wrongly assuming that the group extended 
to all azimuthal angles. 

Figure [3] shows the eccentricity and angular momentum 
distributions of each group compared to the distributions 
estimated in the parent populations. This figure shows that 
the stars in a single group do not have the same angular 
momentum. The rotation period of a star in the galaxy can 
be estimated from a star’s angular momentum. A spread in 
angular momentum in the group implies differential rotation 
between the higher and lower angular momentum members 
of the group. We consider how long it would take a disrupted 
cluster, with stars originally at the same azimuthal angle but 
with different angular momenta to shear out so that stars 
are located at every azimuthal angle in the Galaxy. 

Because the angular rotation rate Q, ~ v^/L is approx¬ 
imately inversely proportional to the angular momentum, 
(with v c the circular velocity and for an approximately flat 
rotation curve), the time it takes an initially compact clus¬ 
ter with a spread in angular momentum values dL to shear 
by 27r in azimuthal angle is 

a '* p ii 

with P the mean rotation period of the cluster. The rota¬ 
tion period at the Solar neighborhood is ~ 0.24 Gyr. In 
4 Gyr there have been approximately 16 rotation periods 
giving dL/L ~ 0.06 for a group that has sheared by 2 tt 
and is now distributed at all azimuthal angles. Using a so¬ 
lar value of Llsr — 1870 for the cluster mean, we estimate 
dL = 120 km s -1 kpc is required for the cluster to shear 
to 27r at At = 4 Gyr. A parent population with a distri¬ 
bution with dispersion dL < 120 km s -1 kpc and age of 
4 Gyr would not be evenly distributed in azimuthal angle. 
However, older populations with larger angular momentum 
dispersions would be evenly distr ibuted in the Galaxy. 

We use the rotation curve by Allen & Santillan ( 199l][ ! 
to compute the azimuthal angle 


5 The second term of equation 5 by |Allen &; Santillan| |l991 ) 
should have the opposite sign. 
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Figure 7. Distribution of parent populations of abundance groups, taking into account the orbital probability and using the eccentricities 
and angular momenta of the stars in the group. Each group can be identified by its age on the top of the plot. Eccentricity and angular 
momentum of stars in each abundance group are plotted as points on top of number density of the estimated parent distribution. Each 
panel shows a different abundance grouping with group properties listed in Table ^ The color bar shows the estimated number of stars 
per eccentricity and angular momentum bin with bin width and height de = 0.01 and dL = 20 km s —1 kpc. The dashed yellow lines 
border the forbidden region. 
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Figure 8. The angular momentum distribution of the |Bensby et a l. ( 201 sample is shown in blue (with axis on the right) and that of 
the Geneva Copenhagen Survey (GCS) in red (with axis on the left) with overlap regions shown in purple. The distributions are shown 
as a function of the V tangential velocity component and using 10 km s —1 kpc bins. The vertical axes show the number of stars in these 
bins. The GCS histogram when multiplied by the function in equation [ 2 ] is shown with black dots (with axis on the right). The black 
dots match the |Bensby et a l. (2014) histogram giving us an estimate of the selection bias compared to the GCS. 



Figure 9. The selection function /(L) (equation [ 2 ]) is shown here as a function of angular momentum in km s 1 kpc. The |Bensby et aT] 
(2014) sample contains proportionally more low angular momentum stars than the Geneva Copenhagen Survey. 


0(L) = Q(L)At (6) 

as a function of angular momentum for a population that is 
initially at the same azimuthal angle at birth. Here At is the 
age of the group and 0 is computed modulo 2ty. After At, the 
more rapidly rotating stars (at lower angular momentum) 
will have increased in 6 more than those rotating slower 
(at higher angular momentum). For three of the groups in 
FigurepT|we show the angle 0(L) as a side panel. For the old 
groups, 0 increases rapidly over a small change in angular 
momentum. As the angular momentum distributions for the 
old groups are large, they are likely to be well distributed in 
the Galaxy. In contrast, the younger groups contain peaks in 
the estimated parent populations that are narrow in angular 
momentum width, and 0 varies relatively slowly across that 
width. In the peaks of the youngest two groups, we may have 
overestimated the parent populations by a factor of a few 
if they are not evenly distributed in the Galaxy. While we 
may have overestimated the number of stars in the youngest 


two groups (and for them only at angular momenta near 
that of the local standard of rest), we have probably not 
overestimated the number of stars in the older groups. 

In this discussion we have neglected phase variations 
in the epicyclic angle. However, the epicyclic frequency is 
faster (about 40% faster) than the angular rotation rate, so 
we expect the shearing in epicyclic angle takes place faster 
than in azimuthal angle. 

Equation [6] assumes that stars were initially at the same 
azimuthal angle and had a similar angular momentum dis¬ 
tribution. Heating and migration could have taken place well 
after the birth of the group. In this case the group would be 
less evenly distributed than estimated using its age and its 
current angular momentum distribution. If the abundance 
group originated in a star cluster that remained bound for 
a long time before disrupting (e.g., 
then the group would be less evenly distributed than esti¬ 
mated here. However as a recently disrupted cluster should 
have a very narrow angular momentum distribution, more 


Lamers & Gieles 2006) 
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Figure 10. Distribution of parent populations of abundance groups. Similar to figure [7] except a bias against high angular momentum 
stars has been removed. Each panel shows a different abundance grouping with group properties listed in Table ^ Eccentricity and 
angular momentum of stars in each abundance group are plotted as points on top of number density of the estimated parent distribution. 
The dashed yellow lines border the forbidden region. On the right for three of the groups, we show an additional panel plotting the 
azimuthal angle 0(L) for a group originally at the same angle but after a time equal to the age of the group. The y-axis on these rightmost 
panels is angular momentum using the same scale as for the abundance distribution and the a>axis is the azimuthal angle 6. Only when 
0(L) is slowly varying should we have overestimated the parent population by assuming an axisymmetric distribution. We have likely 
only overestimated the parent population size for the two youngest groups and only at angular momentum near 1800 km s —1 kpc. 
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recent heating and migration rate would be required to ac¬ 
count for wide current eccentricity and angular momentum 
distributions. 

The parent populations appear to be clumpy, how¬ 
ever this could be due to sparse sampling. Alternately 
phase wrapping due to shearing of azimuthal and epicyclic 
variations could also cause clumping along this boundary 
(Minchev et al. 2009). To estimate the parent population 
distributions we divide by a probability that is sensitive to 
the eccentricity and angular momentum value near the for¬ 
bidden boundary (as we can see from the sampling we used 
in Figure [ 4 ). Along the forbidden region boundary, a small 
error in eccentricity or angular momentum could give a dif¬ 
ference in probability of a factor of a few, and it is precisely 
in this region where most of the stars are located because 
that is the only region where the probability of finding a 
star is high. Errors in eccentricity and angular momentum 
measurements could cause the appearance of clumping near 
the forbidden boundary. We have minimized this effect by 
taking the maximum probability within the estimated er¬ 
rors for each data point. Nevertheless a small variation in a 
star’s eccentricity and angular momentum along this bound¬ 
ary causes a large change in probability and we should be 
careful when interpreting structure in the parent popula¬ 
tions. 

We see from Figures [3j[7| and [To] that the youngest two 
groups have low eccentricity means and dispersions. The es¬ 
timated parent populations are large, greater than a million 
stars, the large size arising because a small fraction of the 
thin disk stars were selected for study by Bensby et al. (2014) 
and our correction for this selection increased the estimated 
number of parent stars. Both groups contain weak tails in 
the distribution extending to higher eccentricity. Only the 
4.8 Gyr old grouping exhibits a tail toward higher angu¬ 
lar momentum, corresponding to stars coming from outside 
the Solar galactocentric radius. It is difficult to determine 
whether the parent population distribution has a large an¬ 
gular momentum dispersion (width) as the mean angular 
momentum values are near that of the Sun and low eccen¬ 
tricity regions above and below this value lie in the forbid¬ 
den region. If there was a large low eccentricity population 
just interior to the Sun, then the higher eccentricity tails 
suggest that the eccentricity width of the parent population 
is wider at lower angular momentum than near L ~ 1800 
km s _ 1 kpc. The estimated parent distributions suggest that 
most stars in the parent populations have not significantly 
migrated (changed in angular momentum) in the last 4-5 
Gyr, though the tails in the parent distributions are sig¬ 
nificant. Perhaps the same population that migrated also 
increased in eccentricity dispersion and a skewed Gaussian 
model for migration might be preferred (see Figure 3 by 
Bland- Hawthorn et al.|20l0 ). 


The 7.1 year old group has a moderate width in its an¬ 
gular momentum distribution with a standard deviation of 
200 km s -1 kpc and a mean of L ^ 1640 km s -1 kpc(and 
for comparison to the other groups see Table |2j). The mean 
angular momentum value is below that of the local stan¬ 
dard of rest, and much of the parent population lies dis¬ 
tant from the forbidden region, though the parent popu¬ 
lation could extend to lower angular momentum (L < 1500 
km s -1 kpc) and low eccentricity (e < 0.2), and into the for¬ 
bidden region. The tail of the distribution below L = 1300 


km s -1 kpc and e > 0.3 suggests that the parent popu¬ 
lation could contain low eccentricity stars below L = 1300 
km s -1 kpc, as the eccentricity dispersion there is larger 
than at the mean L « 1640 km s -1 kpc. The parent angular 
momentum distribution (shown in Figure [ 3 ]) has one strong 
major peak, similar to those of the two youngest groups. 
In contrast the three oldest groups have much wider angu¬ 
lar momentum distributions (also see the standard devia¬ 
tions listed in Table [ 2 ]) . The shape of the parent population 
angular momentum distribution for the 7.1 Gyr old group 
suggests that many stars have not significantly migrated, 
however both width and fraction of stars in the low angu¬ 
lar momentum tail are higher at 7.1 Gyr than for the two 
younger groups. 

The peaks in the parent populations of the three 
youngest groups suggest that the bulk of their stars experi¬ 
enced little migration within 7 Gyrs. Tails in the distribu¬ 
tions imply that stars that have migrated in these groups 
have also increased in eccentricity dispersion. However, a 
thin disk group that increased in angular momentum dis¬ 
persion (due to migration) without increasing in eccentricity 
dispersion would not have stars present in the solar neigh¬ 
borhood unless its mean angular momentum was near that 
of the LSR. 

As none of the peaks in the distributions (Figure |10| ) 
for the 3 older groups contain many stars we don’t attribute 
any significance to the individual peaks. However, the oldest 
groups have both wide eccentricity and angular momentum 
parent distributions, suggesting that both heating and mi¬ 
gration has taken place. 


5 SUMMARY AND DISCUSSION 


We summarize our primary findings here. A discussion fol¬ 
lows. 

1. We find that stars in the 6 largest abundance groups 
by Mitschang et al. (2014) tend to lie near a boundary 


in angular momentum vs eccentricity space where the 
probability is highest for a star to be found in the So¬ 
lar neighborhood, assuming a relaxed parent population 
evenly distributed in azimuthal and epicyclic angles. The 
stars that are most likely found are those with orbital 
apocenter approximately equal to the Sun’s galactocen¬ 
tric radius. The bias has previously been described as a 
crossing time bias ( Mayor et al.||l97 7). 

2. Using the probability for a star to be located in the Solar 
neighborhood (as a function of eccentricity and angular 
momentum) and a crudely estimated selection function 
for the sample, we estimate that the parent populations 
of the abundance groups range from 200,000 to a few 
million members. 

3. The two youngest groups lie nearest forbidden bound¬ 
aries, implying that there could be a significant popula¬ 
tion of group stars that cannot be seen in the Solar neigh¬ 
borhood. However the two youngest groups are the least 
likely to be evenly distributed azimuthally in the Galaxy 
and by assuming an even distribution we may have over 
estimated the size of the parent populations by a factor 
of a few. The angular momentum dispersions of the older 
groups imply that the parent populations are distributed 
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Table 2. Estimated Properties of the Parent Populations 


GID 

Age 

(Gyr) 

Ah 

n 2 

(L) 

(km s - 

°L 

1 kpc) 

<e> 


5 

4.0 

275,550 

3,786,181 

1772 

102 

0.08 

0.06 

1 

4.8 

612,173 

7,573,869 

1754 

167 

0.12 

0.06 

3 

7.1 

466,881 

4,993,675 

1661 

88 

0.11 

0.07 

4 

10.1 

323,715 

871,146 

1511 

266 

0.25 

0.11 

2 

10.2 

315,132 

1,105,883 

1601 

310 

0.23 

0.11 

6 

12.1 

80,620 

267,500 

1679 

282 

0.23 

0.15 


GIC is the group number given by |Mitschang et ah'] ([2014]) and their estimated 
age given in Gyr. N± is the estimated size of the parent population computed 
using equation^ and taking into account the probability of detecting an orbit 
in the solar neighborhood. N 2 is the estimated size of the parent population 
computed using equation and in addition corrects the probability with an 
estimate for the selection function for the observed sample. Mean eccentricity 
and angular momentum standard deviations are computed from the derived 
parent populations (shown in Figure [l0| ). 


at all azimuthal angles in the Galaxy and that we have 
not overestimated the sizes of their parent populations. 

4. Assuming that mean angular momentum is similar to 
that at birth, the width of the parent populations of the 
thin disk groups suggest that the bulk of their stars ex¬ 
perienced little migration within 7 Gyrs. Tails in the dis¬ 
tributions suggest that stars that have migrated in these 
groups have also increased in eccentricity dispersion. In 
contrast, the parent populations of the thick disk groups 
exhibit both wide angular momentum and eccentricity 
distributions suggesting that both heating and radial mi¬ 
gration has taken place. 


Here we assumed that eccentricity and inclination dis¬ 
tributions are not correlated and have ignored the vertical 
motions. Using the vertical velocities it is possible to esti¬ 
mate the inclination distribution of the parent populations. 
Stars with high inclination are less likely to be detected 
within lOOpc of the Sun (Mayor et al. 1977) and we have 
not taken this into account in our estimate of the parent 
populations. The numbers of stars in the older groups, with 
the highest vertical amplitudes, have been underestimated 
by a factor of a few due to this neglect. 

A large cluster may self pollute with supernova and so 
may not remain chemically homogeneous. Consequently, sin¬ 
gle abundance populations ar e estimated to have sizes below 
2 x 10 5 stars (section 3.2 by Bland-Hawthorn et al.pQlO ). 
The large sizes for the parent populations estimated here 
are a concern as they are above this limit. One possibil¬ 
ity is that each group may be comprised of similar but not 
identical fragments (Mitschang et al. 2014). Or the large 
groups may be part of a co-eval population composed of 
stars born nearly at the same time , and with similar abun¬ 
dances (Blanco-Cuaresma et al. 2015), but not necessarily 
all born in the same place (a co-eval but not necessarily 
co-natal parent population). Alternatively the large parent 
population sizes could be attributed to overestimation re¬ 
sulting from our assumption of an axisymmetric and mixed 
parent distribution. Smaller abundance groups were found 
by Mitschang et al. (2014) and these would be consistent 


with the smaller parent sizes estimated for chemically ho¬ 
mogenous populations. 


We mention some uncertainties that affect this study. 


Bensby et al. (2014) gave no selection description for their 


entire sample. We crudely modeled the Bensby et al. (2014) 
sample distribution by comparing it with the GCS, however, 
the GCS sample itself is taken from two different magnitude 


limited source catalogs and is only complete to 40 pc (Nord¬ 


strom et al. 2004). Future attempts to study parent pop¬ 


ulations of abundance groups will be more robust if they 
are based on well characterized samples, and well character¬ 
ized samples would allow more robust estimates of parent 
populations. 


In this study we used a Monte Carlo simulation tech¬ 
nique to estimate the probability that an orbit family would 
be detected in the Solar neighborhood. We then used this 
probability distribution and the stars in each group to esti¬ 
mate the parent population distributions. However, different 
distributions for the parent populations could be assumed 
from the start and Monte Carlo simulations used to predict 
the number and distribution of stars detected in the Solar 
neighborhood. This approach might alleviate some of the 
difficulties caused by the sparse sampling resulting from the 
few stars in each group. 


As did Bensby et al. ( 

2014 

), we adopted the Milky Way 

model by Allen & Santillan ( 

1991). This study could be 


redone with different or updated Milky Way mass distri¬ 
butions to see how the estimated parent populations are 
dependent upon the underlying assumed Galactic mass dis¬ 
tribution. Both accurate space motions and a good Milky 
Way mass model are needed to better estimate the parent 
population distributions, particularly for stars near the for¬ 
bidden boundary where the probability is a strong function 
of eccentricity and angular momentum. 


We assumed a sharp edged spherical boundary at 100 
pc from the Sun for the solar neighborhood sample. However 
approximately 5% of the stars from the 6 abundance groups 
are at larger distances. Future work could study the impact 
of a selection function that depends on distance from the 
Sun. Errors in measurement of eccentricities and angular 
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momenta have been neglected from this study. These too 
could be more accurately modeled. 

Azimuthal structure in the probability distributions has 
been ignored in this study, however the probability distri¬ 
butions could be sensitive to position with respect to the 
Galactic bar, spiral arms and other dynamical structures 
such as the Galactic warp. 


Future studies may detect variations in the orbital prop¬ 
erties of the groups in different directions allowing a study of 
azimuthal variations and correlations between orbital prop¬ 
erties as a function of distance from the Sun. As more stars 
are identified in a single group, it will be possible to de¬ 
termine whether clumps in e, L are real. Clumps along the 
region of high probability in e, L space might arise because 
of a non-uniform distribution in epicyclic amplitude (e.g., 

: Minchev et al.|2009 ). The location of peaks in the distri¬ 


bution might depend on distance from the Sun, particularly 
if the group is not well mixed in the Galaxy. Detected struc¬ 
tures would be exciting to study with models of how groups 
evolve as they move in the Galaxy. 


In summary, we were surprised by the large sizes of our 
estimated parent populations. The large sizes imply that 
large abundance groups found in the vicinity of the Sun are 
unlikely to be co-natal populations unless they are unevenly 
distributed in the Galaxy. If the groups are not co-natal then 
they may not be comprised of stars exactly the same age. 
The color magnitude diagram fits to the abundance groups 


were no worse than those of open clusters (Mitschang et al. 


2014), suggesting that if there is an age spread in each group, 


it is not large. However, an age spread in the stars in the low 
metallicity groups might contribute to the large eccentricity 
and angular momentum dispersions of these groups. Like¬ 
wise the higher metallicity groups may have lower angular 
momentum and eccentricity dispersions simply because they 
are comprised of younger and thin disk stars. Despite these 
concerns, the increasingly large samples of stars with accu¬ 
rate abundance measurements (e.g., De Silva et al. 2015) 
should be used to study groups of stars with similar abun¬ 
dances and may be used to probe mechanisms such as mi¬ 
gration. However, constraints on the dynamical evolution of 
stellar sub-populations will require larger and better char¬ 
acterized samples, samples that extend away from the solar 
neighborhood and comparisons between observed and pre¬ 
dicted distributions of many streams and groups. 
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